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Abstract- In this paper, Laplace decomposition algorithm 
(LDA) is introduced for the approximate solution of the 
system of homogeneous linear partial differential 
equations. The technique is described and illustrated with 
some numerical applications. The results assert that this 
scheme is rapidly convergent and quite accurate by which 
it approximates the solution using only few terms of its 
iteration scheme. 
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1. INTRODUCTION 


System of homogeneous linear partial differential 
equations has attracted much attention in a variety of 
applied science because of their wide applicability. These 
systems were formally derived to describe wave 
propagation to model the shallow water waves[1-3] and 
to examine some chemical reaction —diffusion model of 
Brusselator[4]. 


In this work, we used Laplace decomposition 
algorithm(LDA) because this scheme provides the 
solution in a rapidly convergent series with componenis 
that are elegantly computed. The Laplace decomposition 
algorithm was first proposed by Khuri[5-6] which is 
further used by Yusufoglu[7] to solve Duffing equation 
and Elgazery[8] for Falkner-Skan equation. The 
modification of Laplace decomposition method 
introduced by Hussain and Khan[9].Khan and 
Gondalf 12] applied Laplace decomposition method fora 
new analytical solution of foam drainage equation. 


Khan and Hussain[13] applied Laplace 
decomposition method on semi-infinite domain. The 
restriction and improvements of Laplace decomposition 
method was given by Khan and Gondall{ 14]. Zafar et.al. 
used Laplace decomposition method to solve Burger’s 
equaton[11]. It is worth mentioning that the proposed 
method is an elegant combination of Laplace transform 
and decomposition algorithm. The advantage of this 
proposed method is its capability of combining two 
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powerful methods for obtaining exact solution. The aim 
of this work is to establish exact solutions or 
approximate solutions of high degree of accuracy for the 
system of homogeneous linear partial differential 
equations. 


II. LAPLACE DECOMPOSITION 
ALGORITHM 


In this section, we present Laplace decomposition 
algorithm(LDA )for solving the system of homogeneous 
linear partial differential equations written in an operator 
form 


Du+K, (u,v)=0, dd) 
Dy+K,(u,v)=0 
with initial conditions 
u(x.0)=f(x), 
v(x,0)=9(X), (2) 


where D, is considered a first order partial differential 
operator, K, and K, are linear operators. 
The method consists of first applying the Laplace 
transform to both sides of equations in system ( 1), we have 
L{D,u} +L{K,(u,v)} =0 
L{D,v} +L{K,(uv)} =0 GB) 
Using the differentiation property of Laplace transform 
and initial conditions (2) in(3), we have 
Lu} =fxy/s— V/s LE Kv) 
Liv} =g(x)/s—I/sLt K,(u,v)} 
Operating inverse Laplace transform on both sides of (4). 


(4) 


we have 
u(x,t)= fx) -L"{1/s LE K,(uy)}} 
v(x,t)=g(x)-L"{L/sL{ Kv) }} (5) 
The Laplace decomposition algorithm(LDA) defines 
the solutions u(x.t) and v(x,t) by the infinite series 
u(x,t) =Z%, u, and v(x ,t) = Dv, 
Substituting (6) in (5), we have 
Seou, =f(x) -L'{ Ws LAK (Zeou, 
Leo Vat 
iy V= B(x) —L"{ V/s LEK, (Zhou, 


Zee Va) }} 


(6) 


(7) 


21 


0 


22 


In general, the recursive relation is given by: 
u, = f(x) and v, = g(x) 

U(X) =—L"{ I/s L{ K,(u,.v,)}}. n20 

and (9) 
V,(X,t)=-L"{ I/s Lf K,(u,,v,)}}. 020 


Il, APPLICATIONS 


In this section, some applications are given in order to 
demonstrate the effectiveness of Laplace decomposition 
algorithm (LDA) to solve the system of homogeneous 
linear partial differential equations. 


A, APPLICATION: 1 


Consider the system of homogeneous linear partial 
differential equations 
uty,-(u-v)=0 


Vv, —ut (u-v) =0 (10) 


with initial conditions 
u(x,0)=e", v(x.0)=e" (i) 
Taking the Laplace transform on both sides of (10) and 
using the differentiation property of Laplace transform 
and initial conditions (11), we have 
L{u}=e'/s - /sL{v,} +1/sL{u-v} 
Liv} = e*/s+I/sL{u,} -I/sL{u-v} 
Operating inverse Laplace transform on both sides of 
(12), we have 

u(x,t)=e'-L" {I/sL{v,} } +L" {1/sL{u-v}} 

v(x ,t)=e"+L"{1/s L{u,}} -L"{1/s L{u-v}} (13) 


(2) 


Then, by using the Laplace decomposition 
algorithm(LDA) which defines the solutions u(x st)and 
V(x ,t) by the infinite series as: 


u(x,t) =X, w, and v(x,t) =ZF,v, (14) 
and the terms u, and y, by the infinite series 
u(x,t) = Ee 5u,, and v,(x,t) = pane (15) 
Substituting (14) and (15) in (13), we have 
=e-L" (Ws Lf Ey, +L {l/s 
Tron} } (16) 


+L" {US L{D5 Qu, }}sL" 
{l/s LE Zi gu, -Drov,}} 


From (16) our required recursive relation is given by: 
u,=e* and 
v=e" 


a7) 


Uy(Xt) = -L" {1/s L{v,.} } + L"{1/s L{u,-v, }}, 
n=0 and 
Vooi(Xat) = L' {1/8 Lfu,.} } -L" {V/s Luv, }}, 


n>0 and (18) 


The first few components of u,(x ,t) and v(x.) by using 


recursive relation (18) as follows immediately 


(19) 
(20) 


u,=te’ and v= te” 
u, =e" t/2! and v, = e"t/2! 


and so on for other components. Using (14), the series 


solutions are therefore given by (21) 
u(x,t) =e"(1 + t+0/2! 4 1/3!+......), 

VOX) = e* (1 +t 7/2! $1338! 0...) 

that converges to the exact solutions 

u(x,t)=e"' 

v(x,0 =e""' (22) 
B. APPLICATION:2 

Consider the system of homogeneous linear partial 
differential equations (23) 
U,- V,+ (uty) =0 

v,—u,t (u+v) =0 

with initial conditions 

u(x,0) = sinhx, v(x,0) = coshx (24) 


Taking the Laplace transform on both sides of (23) then, 
by using the differentiation property of Laplace 
transform and initial conditions (24), we have 
L{u} =sinhx/s+ I/sL{v.} -1/sL{utv} 
L{v} = coshx/s+1/sL{u,} -1/sL{ut+v} (25) 
Operating inverse Laplace transform on both sides of 
(25), we have 
u(x,t) = sinhx +L" {l/s L{v,}} - L"{1/s L{u+v}} 
v(x,t) = coshx +L" {1/s L{u,}} - L"{1/s L{utv}} 

(26) 
The Laplace decomposition algorithm(LDA) defines the 
solutions u(x ,t) and v(x ,t) by the infinite series 
u(x,t) = Diu, and v(x,t) =D" ,p, (27) 
and the terms u, and v, by the infinite series 


U(X,t) =Z cu, and v,(x,t) =LTv,, (28) 


Substituting (27) and (28) in (26), we have 
Zio u,= sinhx + L” {1/s L{ D54v,.}} 

-L*{ l/s L{E2, u,+ E2v,}} 

Zyov.=coshx + L” {1/s L{ Zz ot, } 

-L" {ls L{Z7, u,+ Er0v,}} (29) 
From (29) our required recursive relation is given by: 
u, = sinhx and v, = coshx (30) 


u(x) =L' {l/s Ltv.}}-L"{1/s Lfuty, }}. n20 
and 
V,-(X.t) = L“{ I/s Lfu, 


at} L"{M/s L{u,ty, }}, 120 
Bl) 
The first few components of un(x,t) and vn(x,t) by using 


Tecursive relation (31) as follows immediately 


u,=-tcoshx and — v,=-tsinhx (32) 
u.=t/2!sinhx andy, =t/2!coshx (33) 
u,=-t/3!coshx and —-v,=-t'/3!sinhx (34) 


and so on for other components. Using (27), the series 
solutions are therefore given by 

u(x,t) =sinhx(1 +072! + 0/4! + 
coshx(t+t/3!+t/S!+ 
v(x,t)=coshx(1 +t'/2!+t'/4!+. 
sinhx(t+t/3! + /5!+ 
that converges to the exact solutions 

u(x,t) =sinh(x-t) and v(x,t)=cosh(x-t) 


IV, CONCLUSION 

In this paper, we have successfully developed the 
Laplace decomposition algorithm (LDA) for the 
solution of the system of homogeneous linear partial 
differential equations. The given application showed 
that the exact solution have been obtained even with just 
the two (application:1) or three (application:2) first 
terms of the LDA solution, which indicates that the 
proposed method LDA need much less computational 
work. The proposed scheme can be applied for the 
system more than two homogencous linear partial 
differential equations. 
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